Elastic moderation of intrinsically applied tension in lipid membranes 
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Tension in lipid membranes is often controlled externally, by pulling on the boundary of the 
membrane or changing osmotic pressure across a curved membrane. But modifications of the tension 
can also be induced in an internal fashion, for instance as a byproduct of changing a membranes 
electric potential or, as observed experimentally, by activity of membrane proteins. Here we develop 
a theory which demonstrate how such internal contributions to the tension are moderated through 
elastic stretching of the membrane when the membrane is initially in a low tension floppy state. 



PACS numbers: 87.16.dj, 68.03.Cd 



I. INTRODUCTION 

Lipid membranes are ubiquitous in biological systems. 
Their primary function is to separate cells and organelles 
from their surrounding environments, but they also par- 
ticipate actively in many biological processes. For in- 
stance nerve conduction, exocytosis and endocytosis, and 
production of ATP [l[. 

An important mechanical quantity for membranes is 
tension. For instance endocytosis can be regulated by 
changes in tension Q. Often the tension is controlled 
by external means. For example in micropipette experi- 
ments, where a vesicle is aspirated to a pipette, the ten- 
sion is set through Laplace's law by the pressure differ- 
ence between the inside and the outside of the pipette 
[H . Using Langmuir-Blodgett troughs the tension in lipid 
monolayers on a water surface can be lowered by com- 
pression Or red blood cells can burst by increasing 
osmotic pressure when the induced tension exceeds the 
tensile strength. However, there are also intrinsic ways 
for the tension of a membrane to be modified when the 
membrane is attached to a fixed frame (or encloses a fixed 
volume). Different scenarios for this case are: (i) Changes 
in an electric membrane potential, say when an action po- 
tential travels along a nerve, will induce a change in the 
tension 0, Q . (ii) For the force-dipole model introduced 
in [?| to explain micropipette experiments on membranes 
with active ion pumps, it was argued in @ that the model 
implies an additional active contribution to the tension 
of the membrane, (iii) If a membrane becomes charged 
(say due to a change of pH) then this will also lead to a 
contribution to the tension [i| [l(| ■ (iv) If the membrane 
exchanges material with a reservoir, say lysolipids with 
the surrounding bulk fluid fill ], then the tension of the 
membrane is also likely to be affected. However, this sit- 
uation will be more complicated than the previous three 
cases. For instance, the overall equilibration of the sys- 
tem will depend on balance equations for the chemical 
potentials of the lysolipids in and out of solution. We 
will not consider this case further in the following. 

For each of the situations (i)-(iii) theoretical predic- 
tions have been made for the additional contribution that 
the effect will add to the tension. The question we will 



address in this paper is: how will the total tension of 
the membrane respond to such contributions. The ques- 
tion has gained actuality due to recent experiments on 
the fluctuation spectrum of active membranes [l2j. In 
these experiments a decrease in tension from typically 
4 x 1CU 7 N/ni to 0.5 x 10 -7 N/m was observed when ion 
pumps in the membrane were activated. The magnitude 
of this change is much smaller than the a priori contri- 
bution ~10 _3 N/m predicted by the force-dipole model 
of the membrane protein activity @, Q ■ From the vast 
discrepancy between these numbers one might be lead to 
conclude that a more complicated force distribution than 
a dipole must be used to explain the data. However, in 
Q arguments were sketched indicating that the a priori 
contribution would not necessarily be equal to the total 
change in tension. The purpose of the present paper is to 
present systematic arguments following the formalism of 
[T3I . [hi ] to show how elastic contributions will moderate 
the effect of intrinsically applied tension. We will limit 
ourselves to situations in which only the tension is mod- 
ified, and not other parameters such as bending rigidity, 
compressibility etc. 

The paper is organized as follows. We will first present 
our model in Section [TXJ Then we will calculate fluctua- 
tions and the behavior of the tension in Section HID In 
Section [W] we will discuss the consequences for the dif- 
ferent examples of intrinsic tensions mentioned above. 
Finally we give conclusions and outlook in Section [Vj 

II. THE MODEL 

As our starting point we will take a Hamiltonian with 
three contributions 

U = U C + U S + U &AA . (1) 

The first contribution is the Helfrich bending energy (l5| 

U c = J dA^{2Hf (2) 

where the integral is over the area of the membrane 
and H is the mean curvature at a specific point on the 
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membrane. In Monge gauge, where the membrane is 
parametrized by the height z = h(x, y) of the membrane 
above the xy-plane, we have to second order in h: dA~ 
dxdy{l + [{d x h) 2 + (dyh) 2 }/2} and {2H) 2 = (d 2 h + d 2 h) 2 . 
The second contribution represents elastic stretching con- 
tributions. Expanding around a prefered area A$ one can 
write 
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where A = J dA 1 is the actual area of the membrane and 
K a is the area expansion modulus. These first two contri- 
butions have been found to describe well experiments on 
lipid vesicles, including micropipette experiments span- 
ning a wide range of tensions [ill E3 ■ We have not in- 
cluded a spontaneous curvature or related elastic contri- 
butions (as in the ADE model [l6|). since they will not 
matter for the discussion below (in the Monge gauge H 
is a total derivative up to second order in h and thus 
maximally contribute a boundary term). Similarly, the 
Gaussian curvature is not included since by the Gauss- 
Bonnet theorem it also only contributes a boundary term 
when the topology of the membrane is fixed. The final 
contribution in Eq. (TTJ) represents the additional intrin- 
sically applied tension 



The notation introduced here is that for a quantity X 
its thermal average is: (X) = TvX exp(— ^~H')/Z = 
J[dh]Xexp(-(3n')/Z, where Z = J[dh] exp(-0H') is 
the partition function and [dh] a measure for integrating 
over all possible shapes. 



III. FLUCTUATIONS AND TENSION IN 
MONGE GAUGE 

To connect with the tension obtained in for instance 
a video microscopy experiment we need to find the 
magnitude of the fluctuations of the membrane shape. 
To calculate these we choose the Monge gauge intro- 
duced in the last section, and assume that the mem- 
brane is attached on a quadratic frame of side length 
L. Expanding in a Fourier series with coefficients h q = 
J dxdy e~ l( - xq * +yqy ' > h(x, y), where q x = 2nn x /L and q y = 
2-KUy/ L with n x and n y integers, we have to second order 
in h and its derivatives 
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The question we will investigate in the following is: "how 
does the intrinsically applied tension cr a dd affect the fluc- 
tuation spectrum?" 

An obstacle to answering this question is the nonlocal- 
ity of T~L S . To get around it we will follow the approach 
of [Hldii and use the Hubbard- Stratonovich transforma- 
tion to write = k B T is Boltzmann's constant times 
temperature) 
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meaning that we introduce a new auxiliary field, J, and 
replace the Hamiltonian by 



M = H C + a add A 
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We then take the thermodynamic limit and evaluate the 
integral over J using the method of steepest descent, i.e., 
we assume that fluctuations in J are small and can be 
neglected (for a justification see [Oj])- Thus we let J 
take its stationary value obeying dF/dJ = 0, where F = 
— P^ 1 lnTrexp(— f3H') is the Helmholtz free energy. This 
gives 



J 



A Q 



~{{A)-A ), 



(7) 



where q = yq 2 +q 2 . Using the equipartition theorem 

one finds within this Gaussian approximation the fluctu- 
ations in the shape to be 
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Defining the tension a as the tension that can be observed 
in video microscopy experiments of fluctuating vesicles, 
i.e., the coefficient in front of the q 2 term in the denom- 
inator of Eq. pT7|) (as done for instance in the study of 
active membranes in [T2| ) we have o = <r a dd + J- Insert- 
ing in Eq. (jSJ) we can obtain the excess area a of the 
membrane relative to the frame 
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To evaluate the sum over n x and n y we will assume that 
the tension lies within the interval K<?min ^ a ^ K 9maxi 
where <7 m i n = 2ir / L while the upper cut-off g m ax ~ 
2ir/(2d) with 2d being the membrane thickness. The 
lower bound corresponds to the experimental require- 
ment that a tension dominated regime is observable in 
the fluctuation spectrum, while the upper bound will not 
be very far from the tensile strength of the membrane. 
With this assumption we can approximate J2 n n ~ 
L 2 qdq/(2n) and get 
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Since fluctuations at high q values are not affected by 
changes in tension we will simply assume that the cut-off 
(/max is also constant. Thus combining Eqs. (fT2)) and flj|) 
we get 

da — doWd = dJ = — - — da = —a c — , (13) 
Aq a 

where a c = Kajfa^jr- Rearranging this we get a re- 
lationship between changes in the intrinsically applied 
tension er a dd and the observable tension a 

(l + ^) da = da add (14) 

This is the main result of this paper. The formula can 
straightforwardly be integrated to give the additional 
contribution as a function of the observed tension and 
an integration constant which we take to be the tension 
<7 ro f in the absence of the additional contribution 

0-add = cr - cr ref + a c In . (15) 

fref 

If we restrict ourselves to the low tension limit a <C a c 
(but still having a 3> nqf nin ) then we can obtain the 
relation 

a = <x re f exp(cr add /cr c ). (16) 

In this limit the additional contribution will not affect 
the observed tension with its full effect. Instead it will 
be offset almost completely by a change of the membrane 
area A leading to a countering clastic contribution dJ 
to the observed tension. Note that we have used the 
assumption that a 3> nqf nin to derive Eq. P^)l . Thus 
we should not expect Eq. (fT6"|) to hold for arbitrarily 
large negative cr a dd- In particular, one should expect a 
to become negative at sufficiently negative cr ac id- And 
the membrane will then become unstable if the tension 
reaches a < —Kq^ in . 

In the opposite limit a 3> a c one simply has 

da = da add . (17) 

In this limit the excess area is too small for the membrane 
to contract and offset the additional contribution to the 
tension. Thus the additional contribution will have its 
full effect on the observed tension. 



IV. DISCUSSION 

From the results derived in the previous section we 
see that additional contributions to the tension will have 
to be comparable in size to a c to have an observable 
effect on the fluctuation spectrum. If we take typical 
values K a 0.2 N/m, A « L 2 , k « 20fc B T Qj] one 
gets a c « 4 x 10~ 4 N/m. Let us try to compare this 
value with the expected values for the first three intrinsic 
contributions mentioned in the introduction. 



(i) If an electric potential V m is applied across a lipid 
membrane then it was found in 0, Q that a negative 
contribution to the tension proportional to will be 
induced. For potentials of the order 50 mV it was found 
that the contribution would be cr a dd ~ — 10~ 5 N/m fl8| . 
Thus with the above estimate for a c this electric contri- 
bution to the tension should not be easily observable. 

(ii) For the case of active membranes, then we see 
that the prediction of a tension contribution of the order 
~10 -3 N/m from the force-dipole model for a membrane 
in a floppy state by Eq. (|16[) results in a change of tension 
corresponding to a couple of factors of e. Alth oug h this 
estimate agrees well with the observations of [12[, then 
the estimate is highly uncertain. A major effect, which 
is being neglected when applying Eq. ([T5| . is the factor 
two change in "effective temperature" (or more precisely 
a coefficient corresponding to the factor kBT/(8irn) in 
front of the logarithm in Eq. ((H])) measured for this 
type of active membrane [l!| . The effect of the increased 
"effective temperature" on the tension should be signif- 
icant since the implied increase in fluctuations at short 
wavelengths will compete for excess area with the fluctu- 
ations in the long wavelength tension dominated regime 
[20j . The theory developed here do not take into ac- 
count such changes in "effective temperature" . Thus the 
only conclusion that can be drawn for the active mem- 
brane case is that one cannot immediately conclude from 
the large difference between the additional contribution 
of ^10~ 3 N/m estimated in Q and the chan ge i n ob- 
served value of order ~10~ 7 N/m reported in [12( that 
the force-dipole model does not hold. Further investiga- 
tions are necessary before we have an understanding of 
how the activity modifies a basic physical property like 
tension for lipid membranes. 

(iii) Another case is if a membrane becomes charged, 
for instance due to a change in pH when some of the 
lipids are acidic. Electrostatic interactions will then gen- 
erate a tension as calculated in for example. The 
additional contributions to the tension can easily reach 
a magnitude around a c even if just one percent of the 
lipids becomes monovalcntly charged under physiological 
salt conditions. However, such a change in charge content 
will also alter the bending rigidity of the membrane sig- 
nificantly, breaking the assumption of unaltered bending 
rigidity made for the derivations in this paper. 



V. CONCLUSIONS AND OUTLOOK 

The present paper studied changes in the tension ob- 
served from fluctuations when so called intrinsic contri- 
butions were applied. It was found that when the mem- 
brane was in a low tension state with excess area then an 
intrinsically applied tension would be almost fully can- 
celled by an elastic contribution arising from contraction 
of the membrane area. 

A major limitation of the present work is that it does 
not consider situations with changes in bending rigidity, 



4 



area expansion modulus or prcfcred area, not to mention 
situations with non-equilibrium activity leading to addi- 
tional sources of random noise in the system. An exten- 
sion to these situations would have to take into account 
how fluctuations change also at very short wavelengths, 
around the scale determined by the upper cut-off for the 
wavenumbers (<7 max )- 

We also did not discuss effects on the mechanical ten- 
sion of the membrane (the tension obtained by differenti- 
ating the free energy of the membrane with respect to the 
area of the frame on which it is attached) . There is still 
some controversy about the relation between mechanical 



tension and the tension observed from fluctuaions (see 
flil [UI23 ) ■ We avoided this issue here by only consid- 
ering a fixed frame area L 2 . 
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